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Abstract 

The stochastic transport of suspended particles through a periodic pattern of obstacles in mi- 
crofluidic devices is investigated by means of the Fokker-Planck equation. Asymmetric arrays 
of obstacles have been shown to induce the continuous separation of DNA molecules of different 
length. The analysis presented here of the asymptotic distribution of particles in a unit cell of 
these systems shows that separation is only possible in the presence of a driving force with a non- 
vanishing normal component at the surface of the solid obstacles. In addition, vector separation, 
in which different species move, in average, in different directions within the device, is driven by 
differences on the force acting on the various particles and not by differences in the diffusion co- 
efficient. Monte-Carlo simulations performed for different particles and force fields agree with the 
numerical solutions of the Fokker-Planck equation in the periodic system. 
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In recent years different separation strategies have been explored for the development of 
integrated microfiuidic devices that could provide rapid and efficient fractionation of par- 
ticles in solution. Particularly interesting are novel methods that take advantage of the 
intrinsic Brownian motion of small particles suspended in solution, and the dependence of 
the magnitude of the thermal motion on the size of the particles, to fractionate the initial 
suspension as it is transported through the microfiuidic device. Following the pioneering 
work by Austin and Volkmuth 1], several experiments and theoretical investigations have 
shown the possibility to rectify the Brownian motion of suspended macromolecules (DNA 
in particular) by driving them through an asymmetric course of micrometer-scale obstacles 
3, 3, y, 0] (see Fig. n for a schematic representation of the microdevice). These novel meth- 
ods have the additional advantage that separation occurs in a direction different (normal) 
to the average motion of the suspended particles, thus allowing for continuous fractionation, 
in contrast to traditional separation techniques for biological macromolecules such as elec- 
trophoresis ^ . However, the original explanation which assumes differences in the diffusivity 
of the suspended particles as the underlying mechanism leading to the observed separation 



has been recently challenged by experimental observations and theoretical consider- 
ations 8,0]. Here we investigate the transport of suspended particles in spatially periodic 
systems by means of the Fokker-Planck equation, as well as direct Monte-Carlo simulations, 
and discuss the necessary conditions for separation to occur. In particular, we demonstrate 
that, in the absence of a driving force with a normal component to the solid obstacles at 
the fluid-solid interface, the asymptotic concentration of particles is uniform and no lateral 
separation is possible. We also show that in the presence of a normal force at the surface of 
the obstacles, which induces a non-uniform concentration of particles, the vector separation 
of particles, in which different species migrate at different angles within the microdevice, 
is driven by differences on the force acting on the particles and not by differences on their 
mobility or diffusion coefficient. 

The stochastic transport of suspended particles through an array of obstacles in the pres- 
ence of an external force can be described by the Fokker-Planck equation for the probability 
density P(x, t), 



-P(x, t) + V ■ I ^i^(x, t) I - ^ V^P(x, t) = 0, (1) 
where x = (x, y) is the position vector in the 2D microdevice, ks is the Boltzmann constant. 
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FIG. 1: Schematic view of a microfluidic device for the continuous-flow separation of suspended 
particles, equivalent to the one proposed in Ref. ^ | and implemented experimentally in Refs. Q, ^. 
The gray rectangles are solid obstacles. A suspended particle, represented by a small sphere is 
transported through the sieving medium by an external field F(x). We show the unit cell used in 
this work, with dimensions Lx and Ly in x and y directions respectively. We also show, on the 
right, the primitive cell of this 2D periodic system in dotted lines. 

T is the temperature, F(x) is the force acting on the suspended particles, which we have 
assumed to be time-independent, and r] is the friction coefficient |lO|. The diffusion coeffi- 
cient is given by the well known Einstein relation D = kBT/rj, and for a hydrodynamic drag 
(friction) on a spherical particle of radius a suspended in a viscous fluid of viscosity fi and 
moving at low Reynolds numbers we arrive at the Stokes-Einstein relation D = ksT/GTCfia. 

In general we are interested in computing the probability density flux describing the 
migration of the suspended particles relative to the flxed sieve, 

J(x, t) = ^P(x, t) - ^VP(x, t), (2) 

and, in particular, to seek to determine whether the direction in which particles migrate 
is different for different species, leading to vector separation jl^ . llS^ . In a periodic array 
of obstacles, such as the one shown in Fig. ^ it is convenient to introduce the reduced 
probability density and flux jl^ . 

+00 +00 

P(x,t)= J2 J2 P{x + KLx,y + kyLy), (3) 

kx=—ao ky=—oo 
+00 +00 

J(x,t)= ^ ^ J{x + kxLx,y + kyLy), (4) 



where k = {k^, ky) identifies the specific cell within the periodic system and L^, Ly are the 
unit cell dimensions (see Fig. Q). The previous sums are expected to be convergent based on 
the normalization condition for P(x, t). Note that we have assumed an array extending to 
infinity in both x and y directions, or equivalently we are neglecting finite size effects on the 
behavior of the transported species. Then, given the linearity of the Fokker-Planck equation, 
the convergence of the previous sums, and the invariance under discrete translations in the 
periodic system, {xf,yf) = (x + kxL^.y + kyLy), it is clear that P(x, t) is also a solution 
to the Fokker-Planck equation. In particular, P satisfies Eq. with periodic boundary 
(and initial) conditions. It is clear that the suspended particles move between different 
cells but, due to the invariance of the system, it is as if a particle leaving a given cell was 
in fact returned at the equivalent point in the opposite face of the same cell, leading to 
the mentioned periodic boundary conditions in both P(x, t) and J(x, t). In addition, for 
transit times that are long enough so that the system reaches its asymptotic state, we are 
interested in the long time behavior of the system, governed by the steady-state solution of 
the Fokker-Planck equation. Note that only the reduced probability reaches a steady state, 
whereas the global probability density function has no meaningful asymptotic state jl^ . 
Given an initial distribution of particles at the inlet of the microfluidic device, the steady 
state is reached for transit times longer than the characteristic diffusive time required for a 
single particle to explore the primitive unit cell. That is, when the transit time for a particle 
moving through cells, tc ~ NL/u, where u is the average velocity of the particle and L is 
the dimension of the unit cell, is larger than the diffusive time required to explore the cell, 
td ^ L'^/D, where D is the diffusion coefficient. The previous condition is satisfied for small 
global Peclet numbers, Pe = uL/ND ^ 1. In previous experimental works the condition 
of steady state distribution, i.e. small Pe, is typically reached after the suspended particles 
are transported through a relatively small number of cells. In Ref. jsl, for example, the 
authors observe lateral separation for 48.5 kbp DNA molecules with Pe ~ (m ~ l/im; 
D ~ 0.64yum^/s; L ~ 6yum) but no lateral separation for 411 bp DNA molecules with 
Pe ~ 0.5/A^ {u ~ Ifim/s, D ~ 12^m? / s and L ~ 6/im). In both cases, the steady state is 
in fact reached after the DNA molecules are transported through a few cells, compared to 
the total number of cells in the device, A^ ~ 1000. 

In the presence of different species, such as suspended particles or DNA molecules of 
different size, and assuming that the dilute approximation in which particle-particle inter- 
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actions are neglected is valid, we can describe the dynamics of the different particles by 
independently solving the steady state solution of the Fokker-Planck equation with periodic 
boundary conditions for each of the species present in the system. In addition, for imper- 
meable obstacles, the probability density has to satisfy the zero-flux boundary condition at 
the surface of the solid obstacles. The governing equations then read, 

V- = V- |^^P~(x)-^VP°°(x)l =0 (5) 
[ 7] r) J 



J°° n 



= (6) 

as 



where ds is the surface of the obstacles and n is the local normal to the sohd surface. 
Once we obtain the steady state solutions, P°°(x) and J°°(x), we can compute the average 
migration velocity (as well as its orientation), by computing the total probability flux across 
the boundaries of the unit cell, that is, 

r^dyJ^, (7) 

^0 

and the analogous equation for {Uy). 

Let us now consider two cases, classified according to the properties of the force field, that 
are relevant to the separation of particles in microfiuidic devices. Case I: A solenoidal force 
field with no force normal to the solid obstacles. This is probably the most common situation 
in experiments. For example, if one uses electrophoretic driving forces with obstacles that 
are impermeable to the ionic solution then the electric field outside the Debye double-layer 
(which is typically only a few nanometers thick) will be parallel to the solid surfaces as 
the ions carried by the electric field move around the obstacles. This is also the case for 
hydrodynamic convection in pressure driven fiows. In fact, in order to describe the motion of 
an advected particle we can replace F/77 by the local fiuid velocity v, and again the normal 
component vanishes at the surface of the impermeable obstacles. Unfortunately, the steady 
state in these cases is a uniform concentration of particles in the unit cell, independent of the 
details of the forcing field or the physical properties of the suspended species (Note that the 
force field satisfies V • F = 0). This is due to the fact that, since there is no force towards the 
obstacles at the solid surface, there is no concentration gradient developing at the interfaces. 
Therefore, the average velocity will be co-linear with the coarse-grained global force apphed 
on the system and no vector separation is possible, 

(U) = i(-V r rdxdyF{A. (8) 
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FIG. 2: Contour plots for the steady-state reduced probability density function P°°(x) inside a 
unit cell of the microfluidic device shown in Fig. ^ The force is spatially uniform and oriented at 
45°, Fx = Fy. The difference between the two plots is the magnitude of the local Peclet number, 
Pe = Fl/ksT: a) Pe = 1.12 and b) Pe = 5.58. 

If the magnitude of the force divided by the mobility of the particles, F/r], is different 
for different species, then the particles would move at different absolute velocities, but 
they will still move in the same direction (relative orientation within the device). Case 
II: Force penetration into the solid obstacles. In Fig. |21 we show the steady-state solution 
of the probability density for a spatially uniform force oriented at 45° {F^ = Fy). It is 
clear that concentration gradients develop close to the sohd obstacles, which increase with 
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the magnitude of the external force F. It is possible, in fact, to render Eqs. (jS)) and © 
nondimensional using the local Peclet number, Pe = {Ui)/D = Fi/lksT), where ^ is a 
characteristic length scale for the obstacles. Figure |21 presents results for two different Peclet 
numbers, Pe ~ 1 and Pe ~ 5, where we used i = d (see Fig. Q). These Peclet numbers 
are representative of some of the recent experiments on the transport of DNA molecules in 
microdevices. For example, in Ref. 0] the authors report good separation for local Peclet 
numbers in the range 2.1 < Pe < 5.1 and in Ref. the authors observe the best separation 
performance for Pe ~ 3.5. On the other hand, the asymptotic solution in the limit of high 
Peclet numbers is singular, with high concentration gradients present in a thin boundary 
layer of thickness ^ ~ (A;bT)/F^. In this limit, the assumption of point-particles will no 
longer be valid for ^ ^ a. 

The migration velocity will in general depend on the specific force field acting on the 
particles, but we can, however, obtain some interesting results by considering the properties 
of the steady state solution to Eqs. and (jH)). First of all, for a given force field F(x), the 
steady-state probability density function, P°°, is independent of the mobility of the particles 
and, therefore, independent of their diffusivity. As a consequence, if we have species A and 
B with mobilities r]^^ and r/^^ their average migration velocities will be proportional to each 
other, (Ua) = {vb/va)(Vb), and there will be no vector separation in this case. In order to 
test this result, we compare the solution to Eqs (0) and ® with Monte-Carlo simulations for 
the motion of suspended particles. We simulated the motion of two species of particles which 
experience the same force field but have different diffusivities, e.g. the case of equally charged 
particles of different size driven by an external electric field. First, we simulated the motion 
under an electric field that does not penetrate the obstacles and obtained, as expected, a 
uniform distribution of particles. Then, we simulated the motion under the action of a 
spatially uniform force and, in Fig. El we compare the steady-state solution presented in 
Fig. |2l with the simulation results. It is clear that the simulations results agree well with 
the steady state solution to the Fokker-Planck equations with periodic boundary conditions. 
We can also see that the asymptotic distribution of the two different types of particles is the 
same, independent of the value of their mobility. Then, particles with an order of magnitude 
different diffusivities will migrate in the same average direction. Fortunately, in the case of 
DNA molecules, although the electrophoretic velocity is independent of size, the electric force 
acting on DNA molecules of different length is not the same (A ~ 0.3e~/A is the effective 
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FIG. 3: Probability distribution function at the two perpendicular lines, X and Y ^ in the unit cell 
(see Fig. The force is spatially uniform and oriented at 45°, and the Peclet number is Pe = 5.58. 
The solid line corresponds to the solution to Eqs. © and © presented in Fig. The open and 
solid points correspond to Monte-Carlo simulations for two types of particles, which only differ on 
their mobility, and thus their diffusion coefficient, by a factor of 10. 



charge per length of the DNA molecule 



15|). Therefore, it could be possible to achieve 



vector separation by creating an electric field with a non-vanishing normal component at the 
solid surfaces using, for example, permeable obstacles. In the previous analysis, however, 
we neglected the hydrodynamic interactions between the suspended particles and the solid 
obstacles, which would lead to a position-dependent mobility Then, for large enough 
particles compared to the characteristic size of the obstacles, there might be separation 
induced by the difference in permeability of particles of different size 
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We showed that the stochastic transport of particles through periodic patterns of obsta- 
cles can be described by the asymptotic steady-state solution of the Fokker-Planck equation 
applied to a unit cell with periodic boundary conditions and appropriate condition at the 
surface of the solid obstacles. We demonstrated that, in order to obtain wector separation of 
particles, a non-vanishing driving force normal to the obstacles is required at the fluid-solid 
interface. We also showed that the underlying cause driving separation is the difference in 
the force experienced by different species and not the differences in the diffusion coefficient. 
These results could provide guidance for future developments in the design of microfluidic 
devices, such as the use of permeable obstacles in combination with electrophoretic forc- 
ing of the particles. In addition, the numerical solution of the Fokker-Planck equation in 
different geometries could then be used to optimize the design of such microdevices. It is 
interesting to note that an analogous formalism to the one used here has been developed to 
study transport in periodic porous media {discontinuous systems) and a nice and thorough 
discussion can be found in chapter 4 of Ref. There, the authors show that the use of 
such method goes beyond the study of particle migration and also allows the direct com- 
putation of the dispersion coefficient, responsible for the broadening of an injected band of 
suspended particles, an important parameter affecting the resolution in separation devices. 
Finally, we mention that we are currently considering the effect of transient dynamics before 
steady state is reached, which could lead to permanent separation in short devices. 
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